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Abstract 

Time evolution of a circular cosmic string loop is investigated by numerically 
solving the field equations for the scalar and the gauge fields consisting of the 
vortex. It is shown that the result agrees with an analytic estimate based on 
the Nambu-Goto action, which supports its validity in analyzing nonstraight and 
rapidly moving strings. 



The cosmic string scenario is one of the major proposed formation mechanisms 
of density fluctuations in the universe 0. Although it is true that quantum fluc- 
tuations produced during inflation provide realization of primordial density fluc- 
tuations more economically f2|, cosmic string scenario is no less important in the 
present situation of cosmology that the type of the dark matter dominating mass 
density of the universe is yet unknown, because the latter can induce structure 
formation even in the hot-dark-matter-dominated universe in which linear density 
fluctuations on small scales are washed away due to free streaming. Note also that, 
contrary to the usual prejudice, it is easy to have cosmic strings compatible with 
inflation ||, which is anyway indispensable to realize globally homogeneous and 
nearly flat spacetime apart from the origin of fluctuations {| . 

The most remarkable feature of the cosmic string scenario is that it has es- 
sentially only one parameter, namely, the line density of a string, /i. It should 
satisfy G/i ~ 10~ 6 in order to meet various observational requirements of large- 
scale structures || . On the other hand, it also suffers from an upperbound imposed 
by the timing analysis of a millisecond pulsar ||, since the stochastic gravitational 
background radiation emitted by string loops can disturb it [0]. 

In the usual treatment of cosmic strings, their motion is described by the 
Nambu-Goto action, which is proportional to the surface area of the string world 
sheet ||. Using the two dimensional metric tensor on the surface, 

Kb = 9nv Q f b , \i,v = 0,1,2,3, a, b = 0,1, (1) 

where X M (£ a ) is the spacetime trajectory of a string parametrized by a timelike 
parameter (° and a spacelike parameter the action is expressed as 



S = -fJ. / J-dethatdCdC 1 . (2) 



The above expression can be obtained for the Nielsen-Olesen vortex line treating 
it infinitely thin and assuming the Lorentz invariance along the string || . Strictly 
speaking, therefore, the Nambu-Goto action is justified only for static and in- 
finitely straight strings. The equations of motion deduced from the action would 
adequately describe dynamics of long strings which induce structure formation 
because their curvature scale is many orders larger than their thickness or the 
Compton wavelengths of the scalar and the gauge fields consisting of the vortex 
line. However, they may not correctly describe the motion of small loops, especially 
near cusp-forming regions, where velocity of the string reaches that of light. 

Two approaches have been proposed so far to improve description of string 
motion. One is to calculate finite-width corrections to the Nambu-Goto action 



with respect to the ratio of the string width to the curvature radius |L0|. The 
other is to employ a new phenomenological energy-momentum tensor for a string 
loop in which the line density and the string tension are allowed to take different 
values depending on position and time [11]. The latter research claims that these 



two quantities take considerably different values near a cusp forming region and 
that in some cases the line density vanishes at a cusp, which may imply snapping 
of a string loop. If this is the case, previous calculation of stochastic gravitational 
radiation background may have to be reconsidered. 

In the present paper we investigate the validity of the Nambu-Goto action with- 
out resorting to perturbative or phenomenological methods. That is, we reproduce 
a string loop in terms of scalar and gauge fields and trace its motion by numerically 
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solving field equations. We compare a result of numerical simulation with an ana- 
lytic calculation obtained solving the equation of motion given by the Nambu-Goto 
action. As a tractable example for the both analyses, we consider a circular loop. 

In order to generate a vortex solution, we employ the Abelian Higgs model in 
the flat spacetime, whose Lagrangian is given by 

C = —F^ + id^-ieA^id^ + ieA^-V^}, 
Ffj,v = d^A u -d u A^ (3) 

VM = \{\<P\ 2 ~v 2 ) 2 . 

The equations of motion for the complex scalar field 0(x, t) and the gauge field 
^(x, t) read, 

□ + 2ieA»d fl - e 2 A fl A" + iefaA")] + A(|0| 2 - v 2 )(j) = 0, (4) 
dpdvA" - DA,, = ie(<f)*d^ - (jtirf*) + 2e 2 A^*(j), (5) 

respectively. We take A = 0.01, e = 1, and v — 1 and impose the Lorentz gauge 
condition, d^A^ = 0, for numerical calculation. 

In order to set up the initial configuration, using the relaxation method fi"2||, we 
first reproduce the static and infinitely straight Nielsen- Olesen vortex 0, in which 
amplitudes of the scalar field and the gauge field are given as a function of the 
distance from the string core. Then it is bent artificially to make a circular loop. 
To be more specific, we calculate the minimum distance to the trajectory of the 
string core at each point and assign the corresponding field amplitudes just as the 
Nielsen-Olesen vortex, while varying the phase of the scalar field uniformly and the 

direction of the gauge field correspondingly. Since we also set = and A^ = 
initially, where an overdot denotes time derivation, the second time derivative of 
both cj) and A^ is nonvanishing then due to bending. But it is small because we 
take the initial radius of a loop large compared with the core radius of a string. 

Making use of the rotational symmetry, the above field configuration is embed- 
ded in the center of a square of 200 2 meshes perpendicular to the string loop with 
the boundary condition diA^ = and di<p = where di denotes spatial derivatives. 
In fact only a quarter of the square is required for numerical computation due to 
the reflective symmetry. We take the core radius or equal to 10 Ax and the 

initial radius of the loop to be 50Ax with Ax being the mesh size. Then time 
evolution of the loop is traced solving equations. (f|) and (§) in terms of the leap 
frog method with a constant time step At = O.lAx. 

On the other hand, taking a cylindrical coordinate ds 2 = —dt 2 +dr 2 +r 2 d8 2 +dz 2 
and (° = t and C 1 = 9, the Nambu-Goto action (|j) for a circular loop is given by 

S = -27T/i J R(t) sjl- R 2 (t)dt } (6) 

where R(t) is the radius of the loop. Then the equation of motion reads 

R(t)R(t) - R 2 (t) + 1=0, (7) 
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which is easily solved as 



R{t) = R 



cos 



r 



(8) 



where we have assumed the loop is at rest with radius R$ initially at t — 0. The 
speed of the loop motion is given by 



\R(t)\ 



sin 



t 

R 



(9) 



which reaches the speed of light at t = |i?o — 78.5Ax (i? = 50 Ax), when the 
loop shrinks to a point. 

Now we are in a position to display our result of numerical calculation and 
compare it with the above analytic estimate Figure 1 depicts time evolution of 
the potential energy density, V[<f>), of the scalar field up to the stage the loop shrinks 
completely. Due to the initial finite second derivative, the string loop is accelerated 
and its radius approaches zero. Figure 2 shows evolution of scalar field amplitude 
in terms of -F(x, t) = 1 — |0(x, t)\ 2 , which vanishes at the potential minima and is 
equal to unity at the string core. As is seen there the field configuration starts to 
be disturbed around t ~ 60 and after the loop has shrunken to a point, it dissipates 
all the energy into scalar and vector waves without bouncing, which is seen in the 
further simulation. This is quite similar to the fate of a cylindrical domain wall 
consisting of a real scalar field RIB"] . 

Finally we compare time evolution of the trajectory of the string core </>(x, t) = 
with the analytic result ([8p, which is summarized in Table 1. As is seen there 
the result based on the Nambu-Goto action agrees with the numerical calculation 
remarkably well up to the stage that the loop radius becomes as small as the core 
radius and that the speed of string motion is as large as 99% of that of light. 
We can interpret this result as a strong support to conventional analyses based 
on the Nambu-Goto action. We also note that the phenomenon that the string 
loop dissipates all the energy without bouncing has occurred simply because we 
have employed a too idealized configuration of a circular loop. It could be avoided 
if we would consider non-selfintersecting loops. Unfortunately, however, it is a 
formidable task to reproduce such a loop in terms of the scalar and the gauge 
fields themselves. 

In summary, we have calculated time evolution of a circular cosmic string loop 
by directly solving field equations for the scalar and the gauge fields. The result 
agrees well with the analytic estimate based on the Nambu-Goto action, so that 
we may conclude that the action is applicable even for microscopically nonstraight 
strings moving with a relativistic speed. 
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Figure Captions 

Fig. 1 Time evolution of V[(f>] on the 100 2 2-dimensional slice of the string loop 
whose radius is equal to 50Ax initially. The upper right corner corresponds 
to the center of the loop and r depicts the radial direction of cylindrical 
coordinates. These figures shows the epochs t — 0, 30, 60, 70 respectively. 
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Fig. 2 Time evolution of F(x, t) on the same slice for the same string as Fig. 1 
at t = 0,30,60,70,120,160. 

Table 1 





Loop Radius 




Time 


Numerical 


Analytic 


Speed 





50 


50 





10 


49 


49 


0.199 


OA 

20 


44 


44 


r\ o on 

0.389 


30 


41 


41 


0.565 


40 


34 


34 


0.717 


50 


27 


27 


0.841 


60 


18 


18 


0.932 


70 


8 


8.5 


0.985 


75 


2 


3.5 


0.997 



Time evolution of the loop radius in unit of Ax. The third column shows 
analytic estimate based on Eq. (R) and the fourth the speed of the loop motion 

(!)• 
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